Atomic correlations in itinerant ferromagnets: quasi-particle bands of nickel 
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We measure the band structure of nickel along various high-symmetry lines of the bulk Bril- 
louin zone with angle-resolved photoelectron spectroscopy. The Gutzwiller theory for a nine-band 
Hubbard model whose tight-binding parameters are obtained from non-magnetic density-functional 
theory resolves most of the long-standing discrepancies between experiment and theory on nickel. 
Thereby we support the view of itinerant ferromagnetism as induced by atomic correlations. 

PACS numbers: 71.20.Be, 71.10.Fd, 71.10.Ay 



Introduction. Studies of the electronic structure of 
ferromagnetic Ni, in particular by angle-resolved photo- 
electron spectroscopy (ARPES), have since long revealed 
large discrepancies between experiments and results from 
spin-density-functional theory (SDFT) |jl|, |, |, The 
width of the occupied part of the 3d bands is approxi- 
mately Wad = 3.3 eV [||, §, whereas all SDFT calcula- 
tions yield values of M^sdft = 4.5 eV @. In SDFT the 
exchange splitting of majority (t) and minority (|) spin 
bands is large and rather isotropic whereas the ARPES 
results report small and highly anisotropic splittings. 
Typical is the problem of the X21 state which is po- 
sitioned at 0.04 eV below the Fermi energy Ep by the 
ARPES data. In contrast, all SDFT calculations pre- 
dict the X21 state to lie above Ep; so does a study 
using the GW approximation 0|. As a consequence, a 
Fermi surface is predicted with two hole ellipsoids around 
the X point whereas only one has been found in de 
Haas- van Alphen experiments [^j, as later confirmed by 
ARPES |,|,|. 

Here, we present precise ARPES data on the elec- 
tronic structure of Ni, in which full control over the 
three-dimensional wave vector k was achieved using very 
low-energy electron diffraction (VLEED) and ARPES [||. 
For comparison we use the quasi-particle bands from 
Gutzwiller theory jlO[ [ill which is based on a multi-band 
Hubbard model whose electron-transfer amplitudes are 
derived from non-magnetic DFT calculations. We obtain 
good agreement with our ARPES data whereby we re- 
solve for the first time all of the qualitative discrepancies 
between experiment and SDFT for the band structure of 
Ni. We argue that the DFT calculations underestimate 
the ratio between the partial charge densities in the Asp 
and the 3d bands. This appears to be a general problem 
for DFT treatments of transition metals, see Ref. [l^ . 

Experiment. ARPES is a well-established tool for the 
measurement of quasi-particle dispersions £'(k) in crys- 



talline solids. However, a principal problem concerns 
the determination of the surface-perpendicular wave vec- 
tor component k±_ of the photoelectron which is changed 
when the photoelectron exits into the vacuum. Full con- 
trol of the three-dimensional wave vector k within the 
solid can be regained, if the dispersion of the photoelec- 
tron final states is known. The free-electron approxima- 
tion which is usually employed for this purpose [|[ |^, |l^ 
frequently fails because the final-state bands are gener- 
ally far more complicated 

A precise determination of both quasi-particle ener- 
gies and wave vectors is essential for a detailed compar- 
ison between experiment and theory. To this end, we 
have controlled k using the following route. First, angle- 
dependent VLEED is applied to determine the unoccu- 
pied states above the vacuum level whose momentum k 
lies on surface-parallel high-symmetry lines of the Bril- 
louin zone (BZ), relying on the fact that the photoelec- 
tron final states can be viewed as time-reversed LEED 
states. Photoemission via these states, employing the 
Constant-Final-State (CFS) mode, is then used to map 
out the dispersion of the occupied quasi-particle bands 
along these lines. The usefulness and accuracy of this 
method has previously been demonstrated for Cu Q . 

In this work we report on measurements of the Ni(llO) 
surface. The VLEED experiment was performed with a 
conventional LEED unit operated in the retarding field 
mode. ARPES in the CFS mode was carried out at beam- 
line 18-A of the Photon Factory (Tsukuba, Japan); see 
Ref. H for further details. Both VLEED an d ARPES 
measurements were performed along the FX, FY, and 
rS azimuths of the (110) surface, giving access to the 
band dispersions along the TKX, TX, and XL lines of 
the bulk Brillouin zone. The ARPES data on the quasi- 
particle dispersions are presented in Fig. |l|as a gray-scale 
map of the negative second derivative of the photocur- 
rent [|j . Critical point energies are given in Table |. 
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FIG. 1: Grey-scale plot of the negative second derivative of the ARPES intensity for nickel with respect to energy, 
on a logarithmic scale (insets: linear scale) for the TKX and FX directions of the BZ. The dispersionless structure at Ef is 
due to a residual Fermi edge resulting from indirect transitions. Some bands (A2, Xi, E2, E3, Ei towards F) are not seen due 
to unfavorable matrix elements, depending on geometry and chosen final state. Theoretical curves are G-DFT III, see table m 



TABLE I: Binding energies in eV with respect to the Fermi energy Ep (> for occupied states). (...) indicates the spin 
average, error bars in the experiments without spin resolution are given as ±. Theoretical data show the spin average and the 
exchange splittings in square brackets. The second column denote the orbital character of the states, t = t2g, e = Eg, capital 
letters: dominant character. The spin-polarized data (A3;/) were taken at fractions / of the FL distance, with the emphasis on 
the analysis of the exchange splittings. For all calculations, including the SDFT, a spin-only moment of = 0.55/iB was used. 
G-DFT I: A = 12 eV, na = 8.8727; G-DFT II: A = 9eY,na = 8.8727; G-DFT III: A = 9eV, nd = 8.7793. 



Symmetry 


Character 


Experiment Reference G-DFT I 


G-DFT II 


G-DFT III 


SDFT 


(ri> 


S 


8.90±0.30 




8.87 


8.89 


8.86 


8.96[-0.11] 


(r25') 


T 


1.30±0.06 




1.34[0.62] 


1.52[0.62] 


1.44[0.46] 


1.99[0.43] 


(ri2) 


E 


0.48±0.08 




0.62[0.28] 


0.71[0.29] 


0.65[0.195] 


0.86[0.41] 


(^i> 


sE 


3.30±0.20 


i 


2.98[0.57] 


3.39[0.50] 


3.31[0.36] 


4.37[0.20] 


(^3) 


T 


2.63±0.10 




2.60[0.78] 


2.95[0.72] 


2.86[0.54] 


3.82[0.37] 


X21 


E 


0.21±0.03 




0.23 


0.27 


0.165 


0.35 


X21 


E 


0.04±0.03 




0.02 


0.02 


0.01 


-0.09 




T 


0.15±0.03 




0.22 


0.24 


0.10 


0.23 


Ae,(X2) 


E 


0.17±0.05 




0.214 


0.25 


0.155 


0.44 


At.,(X5) 


T 


0.33±0.04 


i 


0.45 


0.49 


0.38 


0.56 


{K2) 


spTe 


2.48±0.06 




2.33[0.70] 


2.64[0.66] 


2.59[0.50] 


3.37[0.32] 




pE 


0.47±0.03 




0.51[0.26] 


0.58[0.28] 


0.51[0.185] 


0.70[0.41] 


{Li) 


sT 


3.66±0.10 




3.21[0.76] 


3.59[0.68] 


3.51[0.515] 


4.56[0.23] 


(i3> 


tE 


1.43±0.07 




1.39[0.48] 


1.58[0.47] 


1.51[0.34] 


2.02[0.40] 




Te 


0.18±0.03 




0.30[0.37] 


0.34[0.41] 


0.215[0.30] 


0.38[0.50] 


{L2'} 


P 


1.00±0.20 




0.32 


0.29[0.0] 


0.97[0.0] 


0.24[-0.12] 


{Wi} 


sE 


0.65±0.10 




0.66[0.28] 


0.76[0.30] 


0.69[0.20] 


0.94[0.39] 




T 


0.15±0.10 


I 


0.23[0.46] 


0.25[0.52] 


0.11[0.38] 


0.23[0.56] 


(A3;l/3) 


ptE 


0.57[0.16±0.02] 


14 


0.64[0.30] 


0.73[0.32] 


0.67[0.22] 


0.90[0.42] 


(A3;l/2) 


ptE 


0.50[0.21±0.02] 


if 


0.54[0.34] 


0.61[0.37] 


0.55[0.26] 


0.76[0.44] 


(^3;2/3) 


pTE 


0.35[0.25±0.02] 


U 


0.35[0.37] 


0.39[0.41] 


0.33[0.29] 


0.49[0.48] 
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Gutzwiller theory. The Gutzwiller theory employs a 
tight-binding Hamihonian for a basis of 4s, 4p and 3c? 
orbitals 



that |<J>) is the ground state of the effective Hamiltonian 
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(1) 



The energy-transfer integrals ti a-;j,(T' between the spin- 
orbitals a at site i and cr' at site j range up to third 
nearest neighbors. They result from a fit to DFT energy 
bands of non-magnetic Ni, calculated with the LAPW- 
WIEN code and using the local-density approxima- 
tion; for details of the fit procedure, see Ref. A value 
of 15meV for the root-mean-square deviation was ob- 
tained for all bands up to 2 eV above the Fermi energy. 
We note that the DFT energy of the 4p-type L2' state is 
0.120 eV below Ep, and the bare 3d band width in our 
tight-binding model is Wbaro ~ W^dft = 4.5 eV. 

The single-particle term is supplemented by the atomic 
Coulomb interactions among the 3d electrons. 
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(3) 



Here, fhi-;r is the projector onto the atomic eigenstate 
\T)i with energy Er- In the spherical approximation, the 
Coulomb parameters U'^^''^^'^'^^''^^ can be obtained from 
three Racah parameters. In this work we choose B = 
0.09 eV and C = 0.40 eV, guided by spectroscopic data 
for isolated Ni^+ and Ni^+ ions; see Table 5.1 of Sugano 
et al. This choice fulfills the 'canonical' ratio C / B k 
4.5, found for transition metal ions both in atomic and 
impurity spectroscopy The Racah parameter A — 

0(10 eV) is adjusted; see below. 

The sizable atomic Coulomb interactions suppress lo- 
cal charge fluctuations which occur in the limit of in- 
dependent particles. Therefore, Gutzwiller ||l^ proposed 
to address the states 



^Gl<f> 



Pn 



(4) 



as a variational approximation for the true ground state 
of H. The real numbers 7?r parameterize the Gutzwiller 
correlator Pq which reduces energetically unfavorable 
atomic configurations in the normalized one-particle 
product state |$). For an open d-shell system, there are 
2^", i.e., of the order of 10^ variational parameters. 

The variational ground state is determined by a mini- 
mization of the expectation value {H)^s,^ with respect to 
the 'internal' parameters rjr and the wave function |$). 
In infinite dimensions ^ this leads to the condition 



(5) 



Sa,a'0<) = %/9^%/g^e<T,o-'(k) + (5a-.<j'(eo. + Ag-) , (6) 

where e^^a'^^) is the Fourier transform of the bare 
electron-transfer integrals U^^-j^^i . The spin-orbital en- 
ergy shifts Act act as 'external' variational parameters 
since they determine |$) via the eigenvalue equation 
{H°^ -E)\<S?) =0. 

The quasi-particle excitations from energy bands given 
by the eigenvalues E(k, ^) oi H"^ Note, that the ef- 
fective single-particle Hamiltonian differs from Hi by the 
renormalization factors < g^- < 1 and the energy shifts 
Act. Therefore, the bare bands of Hi become narrowed 
and mixed in the quasi-particle excitations of the corre- 
lated system. The above results are strictly valid only 
for D ~ 00. However, 1/D corrections were found to be 
small for = 3 single-band Fermi- liquid systems ||ll| . 

The Gutzwiller method overestimates the optimum 
magnetic moment by about 10% to 15%; for a remedy of 
this problem, see Ref. |2^. Therefore, we work with a fixed 
magnetic moment = 0.55/LtB, the experimental spin- 
only moment. Spin-orbit coupling has not yet been incor- 
porated into our Gutzwiller codes as this requires a major 
reprogramming effort. We keep the partial densities Ug, 
Up, and Ud fixed to the values obtained from the DFT 
one-particle Hamiltonian. We have studied other vari- 
ants of applying the Gutzwiller-DFT which allow charge 
flow between Asp and 3d channels, yet the effects on the 
quasi-particle bands are found to be small pO| . 

Results for nickel. In cubic iron-group metals and 
under the limitations of keeping n^, n^, and rip fixed, 
symmetry allows for three independent energy shifts Aa- 
for the 3d electrons. They govern the exchange split- 
tings A{t2g) and A(eg), and the crystal- field sphtting 
^{eg/t2g). For Ni, the exchange splittings of the 4s 
and Ap orbitals turn out to be of very minor importance. 
Therefore, we minimize the variational ground-state en- 
ergy for a fixed magnetic moment /x with respect to three 
external and about 400 internal variational parameters. 

We find the best agreement with the occupied 3d 
band width W^d = 3.3 eV at Xi for a value A = 9eV, 
see table |, column G-DFT II. Then, the condensa- 
tion energy, i.e., the energy gain for the ferromagnetic 
phase as compared to the paramagnetic one, is i^cond ~ 
40meV = O.S/cbIc (Tq = 630 K is the Curie temper- 
ature). This value of A is much larger than in other 
multi-band studies 
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In |gl|, only 3d bands were 
incorporated. If we eliminate in our calculations the 
rather large 4sp-3d hybridizations we also find that val- 
ues oi A — 2 eV . . . 3 eV give ferromagnetic solutions with 
i^cond = 40meV and Wgd « 3eV. The resuhs of [|2), 
where the Asp bands have also been included, are at vari- 
ance with our findings. 
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For all three G-DFT calculations of table |, the en- 
ergy of the state X21 is pinned slightly below the Fermi 
energy, in agreement with experiment. This is an essen- 
tial improvement over SDFT and even SDFT-GW cal- 
culations which predict a second hole pocket at the 
X point. Even a drastic change in the Racah parame- 
ter A from All — Am = 9 eV to — 12 eV does not 
change much the basic features of the Gutzwiller quasi- 
particle band structure. The pinning of X2[ is caused 
by the small variational parameter A(eg), which in turn 
leads to a large fraction of t2g holes in the minority spin 
bands, produced by a relatively large A{t2g)- Because of 
the large nearest-neighbor hopping between t2g orbitals, 
it is energetically favorable to generate as many t2g holes 
in the minority spin bands as possible. 

A more detailed inspection of the table | indicates sev- 
eral significant discrepancies between experiment and the 
results from G-DFT II. Most notably, the energy of the 
state L2', £^011(^2') = — 0.3eV, is much higher than 
in experiment, E{L2') = — l.OeV. This discrepancy is 
caused by the underlying DFT calculation (see also ta- 
ble |, column SDFT), as the 4p-type L2' state is basically 
not affected by correlations in the Sd shell. When we 
lower the bare 4p level in (|^) by 0.75 eV, we consider- 
ably improve the agreement with experiment especially 
for all states near the Fermi energy and thus also for the 
values of the exchange splittings. The photoelectron data 
reveal Ag (-^2) = 170 ± 50meV for the pure d{eg) state 
X2, and AejXs) = 330± 40meV § for the pure d{t2g) 
state X5, which agrees well with the data from the G- 
DFT III. Lowering the Ap level leads to a charge flow of 
about 0.1 electrons from Zd to 4p whereby the number 
of holes in the 'id bands is significantly enhanced. Thus, 
there are more id holes available to generate the mag- 
netic moment, and the exchange splittings are reduced 
accordingly. Note that the G-DFT exchange splittings 



Therefore, we corroborate early ideas by van Vleck p8[ , 
Gutzwiller pof and others on the origin of itinerant fer- 
romagnetism. In contrast to band aspects, they empha- 
sized the importance of local moments in narrow-band 
metals as a consequence of strong atomic correlations. 

Further predictions of our theory, e.g., the increase of 
the exchange splittings as a function of energy, should be 
tested by more detailed spin-resolved experiments. 

W.W. thanks S.G. Louie and M.L. Cohen for their hos- 
pitality during his sabbatical stay. We thank H. Starn- 
berg and P. Blaha for their help with the VLEED data 
acquisition and analysis. This work was supported by the 
Deutsche Forschungsgemeinschaft (446JAP/113/104/0, 
BU1309/2-1, CL124/5-1, WE1412/8-1). 



incr t;ase as a fimctiou uf Lliu binding uuuig^' h}j abouL 
50 %), as can be seen, e.g., from the values of the pure t^g 
states X5, r25', and X^. In contrast, they decrease in 
the SDFT because of the bigger orbital basis used there. 

Conclusions. Hartree-Fock theory was the starting 
point for the Stoner-Slater band theories of itinerant fer- 
romagnetism |4j. For decades these theories have 
been the only ones to provide detailed comparison with 
experiments. In this paper we have shown for the proto- 
typical ferromagnet Ni that, starting from a DFT-based 
one-particle Hamiltonian, the Gutzwiller theory can re- 
solve the main discrepancies for Ni. Wc thereby confirm 
the qualitative findings of previous model studies [|lO[ |2^ 
on the basic mechanism of itinerant ferromagnetism. In 
Ni the id^ , id^, and 3(i^° niultiplets are predominantly 
occupied to accommodate the 3c? electrons' itinerancy. 
However, the occupation of the low-lying id^ spin triplets 
is significantly enhanced over the occupation of the id^ 
singlets. As in magnetic insulators, Hund's first rule is 
found to be an important agent in itinerant ferromagnets. 
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